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ABSTRACT 

The Large and Small Magellanic Clouds (LMC and SMC, respectively) are observed to have charac¬ 
teristic dust extinction curves that are quite different from those of the Galaxy (e.g., strength of the 
2175 A bump). Although the dust composition and size distribution of the Magellanic Clouds (MCs) 
that can self-consistently explain their observed extinction curves have been already proposed, it remain 
unclear whether and how the required dust properties can be achieved in the formation histories of the 
MCs. We therefore investigate the time evolution of the dust properties of the MCs and thereby derive 
their extinction curves using one-zone chemical evolution models with formation and evolution of small 
and large silicate and carbonaceous dust grains and dusty winds associated with starburst events. We 
find that the observed SMC extinction curve without a conspicuous 2175 A bump can be reproduced 
well by our SMC model, if the small carbon grains can be selectively lost through the dust wind during 
the latest starburst about 0.2 Gyr ago. We also find that the LMC extinction curve with a weak 2175 
A bump can be reproduced by our LMC model with less efficient removal of dust through dust wind. 
We discuss possible physical reasons for different dust wind efficiencies between silicate and graphite and 
among galaxies. 

Subject headings: Magellanic Clouds - galaxies: evolution - galaxies: ISM - ISM: dust, extinction 


1. INTRODUCTION 

The Large and Small Magellanic Clouds (LMC and the 
SMC respectively), which are the closest pair of interact¬ 
ing dwarf galaxies around the Milky Way (MW), have 
long served as an ideal laboratory to study many aspects 
of galaxy formation and evolution, such as significant in¬ 
fluences of galaxy interaction on galactic star formation 
and chemical enrichment histories (e.g., Westerlund 1997). 
Structures and kinematics (e.g.. Cole et al. 2005; Nidever 
et al. 2012), age and metallicity distributions of stellar 
populations (e.g. Da Costa & Hatzidimitriou 1998; Piatti 
et al. 2001; Carrera et al. 2008), proper motions (e.g., 
Kallivayalil et al. 2006; Cioni et al. 2014), and gas and 
dust contents (e.g., Kim et al. 1999; Meixner et al. 2011; 
Galliano et al. 2011) have been extensively investigated 
so far. These observed physical properties of the Magel¬ 
lanic Clouds (MCs) have been used to constrain theoretical 
models for the MC evolution, such as the 3D orbits around 
the MW and the LMC-SMC interaction histories and their 
influences on the star formation histories (e.g., Gardiner 
& Noguchi 1995; Diaz & Bekki 2011). 

One of the long-standing problems related to the phys¬ 
ical properties of the MCs is the origin of their extinc¬ 
tion curves that are quite different from those of the MW 
(e.g., Fitzpatrick 1989). The MW extinction curve clearly 
shows the 2175 A feature (“bump”) whereas the 2175 A 
bump is almost absent in the SMC and rather weak in 


the LMC (e.g.. Savage & Mathis 1979; Rocca-Volmerange 
et al. 1981; Prevot et al. 1984; Clayton & Martin 1985; 
Fitzpatrick 1985; Gordon & Clayton 1998; Gordon et al. 
2003). Furthermore, the extinction curves of the MCs 
at shorter wavelengths (1/A > 5/im“^) are different from 
those of the MW in the sense that they more steeply rise 
(e.g., Hutchings 1982; Nandy et al. 1982; Bromage & 
Nandy 1983). These observations suggest that the optical 
and physical properties of dust grains (e.g., size distribu¬ 
tions) in the MCs are quite different from those of the MW 
for some physical reason. 

Pei (1992) adopted the graphite-silicate grain model for 
the MCs and demonstrated that the observed extinction 
curves can be reproduced reasonably well by the mod¬ 
els with different abundance ratios of graphite to silicate. 
Weigartmer & Draine (2001) successfully reproduced the 
observed extinction curves of the MCs and the MW by 
adopting a simple functional form for the size distribu¬ 
tions of carbonaceous and silicate grains. Although these 
previous studies clarified dust compositions and size dis¬ 
tributions required for explaining the observed extinction 
curves of the MCs, it remains unclear whether and how 
the required dust properties can be achieved in the for¬ 
mation histories of the MCs. Since the time evolution of 
dust properties depends on star formation and chemical 
enrichment histories of galaxies (e.g., Dwek 1998), a dust 
evolution model that can explain self-consistently the ob¬ 
served SFH and chemical abundances of the MCs needs to 
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be constructed for better understanding the origin of the 
extinction curves of the MCs. 

Asano et al. (2013) have first incorporated the time 
evolution of dust size distributions into one-zone chemi¬ 
cal evolution models in a fully self-consistent manner and 
thereby discussed what determines the dust size evolution 
in galaxies. Nozawa et al. (2015) have additionally in¬ 
corporated the molecular cloud phase hosting quick dust 
growth and thereby shown that the observed peculiar ex¬ 
tinction curves of high-redshift quasars {z > 4) can be 
reproduced, if the quasars can have both a rather high 
molecular gas fraction (> 0.5) and amorphous carbon in¬ 
stead of graphite. These studies, however, have not at¬ 
tempted to reproduce the observed characteristic extinc¬ 
tion curves with the nearly absent 2175 A feature in the 
SMC and the weak one in the LMC. Although recent nu¬ 
merical simulations of the MCs have shown the time evo¬ 
lution and the spatial distributions of dust abundances in 
the MCs (Yozin & Bekki 2014), they did not show the size 
distribution of dust grains in the MCs. Thus, it remains 
theoretically unclear in what physical conditions the ob¬ 
served extinction curves of the MCs can be reproduced in 
the evolution models of the MCs. 

The purpose of this paper is to investigate the dust ex¬ 
tinction curves of the MCs by using new one-zone chem¬ 
ical evolution models with the evolution of dust abun¬ 
dances, compositions, and size distributions. In partic¬ 
ular, we investigate in what physical conditions the ob¬ 
served characteristic extinction curves can be reproduced 
in the new models. Since star formation histories (SFHs), 
age-metallicity relations (AMRs), and chemical proper¬ 
ties (e.g., [Fe/H] and [a/Fe]) of the MCs have been de¬ 
rived (e.g., Harris & Zaritsky 2009; Rubele et al. 2012), 
we need to construct chemical evolution models that are 
consistent with these observations. Furthermore, recent 
observational studies by infrared space telescopes, Akari, 
Spitzer, and Herschel, have revealed a number of intrigu¬ 
ing dust properties in the MCs (e.g., Meixner et al. 2010). 
These observations can be used to constrain the present 
dust evolution models for the MCs (e.g., Galliano et al. 
2011 ). 

The plan of the paper is as follows: In the next section, 
we describe our new one-zone chemical evolution models 
which include not only some key dust-related physical pro¬ 
cesses (e.g., growth and destruction) but also radiation- 
driven dust wind depending on dust grains. In §3, we 
present the results on the time evolution of the extinction 
curves for the representative SMC, LMC, and MW mod¬ 
els. In §4, we discuss (i) the validity of the present new 
scenario that explains the origin of the observed extinction 
curves of the MCs and (ii) alternative ones. We do not dis¬ 
cuss the details of the chemical evolution of the SMC and 
the LMC and their dependences on the evolution histories 
of the MCs in the present paper, because they have been 
already discussed in our previous papers (e.g., Tsujimoto 
& Bekki 2009; Bekki & Tsujimoto 2010, 2012, BT12). 

Although the formation process of dust in supernovae 
(SNe) and Asymptotic Giant Branch (AGB) stars in the 
MCs have been investigated both observationally and the¬ 
oretically (e.g., Matsuura et al. 2011; Zhukovska & Hen¬ 
ning 2013; Schneider et al. 2014), we do not discuss these 
issues in the present paper. The observed spatial distri¬ 


butions of dust properties in the MCs (e.g., Galliano et 
al. 2011) are not discussed either in the present study, 
because they are beyond the scope of this paper. The 
present one-zone models do not allow us to investigate the 
observed spatial variations in the extinction curves of the 
MCs (e.g., Gordon & Clayton 1998), and we will discuss 
this important issue in our future papers in which hydro- 
dynamical simulations of the MCs with a model for dust 
size distributions will allow for such an investigation of the 
spatial variation of the extinction curves. 

2. THE MODEL 
2.1. Outline 

We adopt the following two steps to derive the dust ex¬ 
tinction curves of the MCs. First, we investigate the time 
evolution of dust grains by using one-zone chemical evo¬ 
lution models of the MCs that have been developed by 
BT12. Second, we derive the extinction curves based on 
the abundances of dust grains by using the “two-size ap¬ 
proximation model” developed by Hirashita (2015; H15). 
In this two-size approximation, the whole grain popula¬ 
tion is represented by small grains (radius a less than 0.03 
/im) and large grains (a > 0.03 /rm), because the major 
grain processing mechanisms work differently on small and 
large grains (H15). The original purpose of this two-size 
approximation method is to use it in statistical (theoreti¬ 
cal) predictions and Nbody-|-hydrodynamical simulations 
of galaxy formation and evolution and thereby to dramat¬ 
ically reduce the total amount of calculation time for esti¬ 
mating the dust extinction curve for each galaxy and each 
particle of a simulated galaxy. We use the method in the 
present study, because we need to run numerous models 
and the method enables us to complete such calculations 
of extinction curves much faster. 

In the first step, each element (e.g., C, O, and Si) is di¬ 
vided into gas-phase metals and dust, and dust is further 
divided into “small” and “large” grains. We accordingly 
investigate the time evolution of the three components, 
gas-phase metals and small and large grains, separately 
for each element (e.g., C, O, and Si) in a model. The 
time evolution of small and large grains can be used for 
the derivation of the extinction curve in the second step. 
As demonstrated by HI5, the two-size approximation for 
a dust size distribution is very useful and convenient, be¬ 
cause a full calculation of the dust size distribution in a 
galaxy is complicated and time-consuming (e.g., Asano et 
al. 2013). Accordingly we adopt the revised version of 
the two-size approximation model (HI5) in order to de¬ 
rive the extinction curve of a galaxy in the second step. In 
the revised version of the two-size approximation model, 
the dust in the ISM of a galaxy is divided into four grains; 
small and large silicate (referred to as small and large “Si- 
grains”, respectively) and small and large carbonaceous 
grains (small and large “C-grains”, respectively). This 
two-size, two-component dust grain model enables us to 
derive the extinction curve of a galaxy based on the chem¬ 
ical evolution of the galaxy, as described later. 

Star formation histories (SFHs) of galaxies can control 
their chemical evolution histories, because chemical evolu¬ 
tion depends on the time evolution of metals ejected from 
SNe and AGB stars, the formation rates of which depend 
on star formation rates (SFRs). Therefore, the dust evo- 
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lution histories of galaxies, which depend on their chem¬ 
ical evolution histories, can be controlled by their SFHs. 
Since we mainly discuss the extinction curves of the MCs, 
we need to choose carefully the model parameters for star 
formation of the MCs. Many observations suggest that 
the MCs have experienced at least a few starburst (SB) 
events in the past 6 Gyrs (e.g., Harris & Zaritsky 2004, 
2009; Rubele et al. 2012, 2015). Furthermore, previous 
and recent numerical simulations of the MCs have demon¬ 
strated that highly enhanced star formation is inevitable 
during the past LMC-SMC-Galaxy tidal interaction (e.g., 
Yoshizawa & Noguchi 2003; Bekki & Chiba 2005; Yozin 
& Bekki 2014). We consider these observational and the¬ 
oretical results in constructing the SFHs of the MCs. 

Since the importance and possibility of dust removal 
from galaxies through energetic radiation-driven dust wind 
in one-zone chemical evolution models have been discussed 
extensively by Bekki & Tsujimoto (2014, BT14), we do 
not discuss this issue in detail in this paper. The present 
model is more sophisticated than BT14 in the sense that 
it investigates the evolution of gas-phase metals and the 
formation and evolution of different dust grains in a self- 
consistent manner. Furthermore, the present model newly 
includes, the accretion of gas-phase metals onto already 
existing dust grains, its dependence on cold gas fraction 
in the ISM, grain growth by coagulation, and grain disrup¬ 
tion by shattering. Thanks to this, we can provide more 
reliable predictions on the time evolution of dust grains in 
a galaxy. 


2.2. Basic equations 


We adopt a revised version of our one-zone chemical evo¬ 
lution models used in BT12 in order to investigate the time 
evolution of gas, metals, and dust in a fully self-consistent 
manner. Although there are only two basic equations in 
BT12, there are four in the present study, which reflects 
the fact that the revised version includes the evolution of 
small and large dust grains. We investigate the time evo¬ 
lution of the total masses of gas (Mg) and that of the total 
masses of metal (M^) and dust (Md) for each element (e.g., 
C and O) in a model. Since the total metal mass includes 
gas-phase metals and dust, the total mass of gas-phase 
metals that is not depleted onto dust grains (M^^g) can be 
estimated as — Md in the present study. The total gas 
mass (Mg^t) accreted onto a galaxy is set to be 1.0 (i.e., 
normalized) in all models of the present study. The ba¬ 
sic equations for the time evolution of gas and metals in 
the adopted one-zone chemical evolution models with dust 
physics are described as follows: 


dMg 

dt 


—ai'ip{t) -b A{t) — w{t) — Wd{t) 


( 1 ) 


dt 


= -aiZi{t)ip{t) -b ZA,i{t)A{t) -b yu,iip{t) 


+yia,i / ^(t - tla)g{tli,)dtla 

Jo 


4” / yagb,i(^agb)^(t tagb)^(tagb)dtagb 
Jo 


-W,{t)-WdAt) , (2) 

where ai is the mass fraction locked up in dead stellar 
remnants and long-lived stars, 'ijj{t) is the SFR, Zi is the 


metallicity for tth element, j/nj, yia.i, and yagb,i are the 
chemical yields for the tth element from type H super¬ 
novae (SNe H), from Type la supernovae (SNe la), and 
from AGB stars, respectively, ZA,i is the abundance of 
heavy elements contained in the infalling gas, and Wi is 
the wind rate for each element. The quantities Ha and Ggb 
represent the time delay between star formation and SN 
la explosion and that between star formation and the on¬ 
set of AGB phase, respectively. The terms (/(Ha) and /lagb 
are the distribution functions of SNe la and AGB stars, 
respectively, and the details of which are described later 
in this section. The term ft-agb controls how much AGB 
ejecta can be returned into the ISM per unit mass for a 
given time in equation (2). The total gas masses ejected 
from AGB stars depend on the original masses of the AGB 
stars (e.g., Weidemann 2000). Therefore, this term hagb 
depends on the adopted IMF and the relation between the 
mass and the lifetime in stars (later described). The term 
Wd,i is the dust wind rate for each element (i.e., Wd is the 
sum of Wd,i for all elements). 

The basic equations for the time evolution of the to¬ 
tal masses of small (Md,s) and large (Mdj) grains in the 
models are described as follows: 


dMd,s,i 

dt 


= -a\DsA'^{t) -b I?s,j(0)A(t) 

^ Md,i,i _ Md,s,i Md,s,^ _ 

'Tsh.,i '^co^i 'Tacc,s,i 


_ Md,s,i 
TSN,s,i 

Wd,s,iit) 


(3) 


dMdj,t 

dt 


-aiDi^^tp{t) -b Ad(0)A(t) -b fc,iEz,i 


Md,i,i 

TSN,l.i 


Mdj,i ^ Md,s,i Mdj,i 

'Tsh,i '^co,i 'Tacc,\,i 


Wd,iAt) 


(4) 


where Md,s,i and Mdy,i are the masses of small and large 
grains for Hth element, respectively, and D\ i are the 
dust-to-gas-ratios for small and large grains, respectively, 
A.i(O) and A,i(0) are the initial dust-to-gas-ratios for 
small and large grains (in infalling gas), respectively, /c,i 
and Eza are the condensation efficiency of dust and the 
ejection rate of a metal element, respectively, tsn.s.i and 
TsN,i,i are the timescales of dust destruction by SNe for 
small and large grains, respectively, Tsh,i is the timescale 
for shattering of large grains, Tco,i is the timescale for coag¬ 
ulation of small grains, and Wd,s,i and Wd,i,i are the dust 
wind rates for small and large grains, respectively. The 
dust wind term Wd,i in equation (2) is therefore equiva¬ 
lent to Wd.si + IFd.i.i- 

The ejection rate of metals from SNe and AGB stars 
(A,i) is described as follows: 

A.I = yii.iV’(f) + 2/Ia,i / ^it - tla)g{tls,)dtu 

Jo 


4” / yagb,i(^agb)^(t tagb)^(tagb)dtagb- (5) 

JO 

The formation rate of dust grains from SNe and AGB 
ejecta at each time step is estimated from and /c,i. 
It is assumed here that the dust produced by stars is in 
the large-grain regime (HI5). It should be also stressed 
that the dust condensation efficiency for each Ah element 
is different between SN la, SNH, and AGB ejecta and 
this difference is properly included in the present one-zone 
models: the way to describe the dust production term in 
equation (4) is not very accurate, and such description is 
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adopted only for convenience. In the implementation of 
the dust production process, the dust production term is 
calculated in our code as follows: 

fc,iEz,i = /c,SNIa,ii^z,SNIa,i + /c,SNII,i^^z,SNII,i 

+ /c,AGB,ii?z,AGB,i, ( 6 ) 

where /c,SNia,i /c,SNii,i and /c,agb,i are the dust condensa¬ 
tion efficiencies for ith element and i?z,SNia,i Ii'z,SNii,i and 
i'z.AGB,! are the metal ejection rates for zth element. 

It should be noted here that the way to describe the 
evolution of gas, metals, and dust in the above equations 
is slightly different from those adopted in BT12 and H15. 
For example, we investigate Mg and instead of the 
gas mass fraction (/g in BT12) and metallicity {Zi) in the 
present study. Also, the timescales for dust accretion, co¬ 
agulation, and shattering are defined for each element in 
the present study: they are assumed to be constant for all 
elements in HI 5. The dust wind for small and large grains 
is separately considered, for the first time, in the present 
study. The details of the models for star formation, chem¬ 
ical evolution, and dust physics in the ISM and the model 
parameters are described in the following subsections. 


2.3. Star formation and gas accretion 

The star formation rate 'ip{t) is assumed to be propor¬ 
tional to the gas mass with a constant star formation co¬ 
efficient and thus is described as follows: 

m=CsfMg{t). (7) 

The normalization factor Cgf is chosen appropriately so 
that the final gas mass fraction (and metallicity) can be 
consistent with the observed one. We assume that Cgf is 
different between (i) the “quiescent phase”, when a galaxy 
shows an almost steady star formation, and (ii) the “star 
burst (SB) phase”, when the galaxy experiences a bursty 
star formation. Many observational studies showed that 
the LMC and the SMC have experienced a few SB events 
in the past 6 Gyr (e.g., Harris & Zaritsky 2006, 2009; 
Rubele et al. 2012). Guided by these observations, we 
investigate models in which a starburst can occur three 
times (“SB models”), and the first, second, third star- 
bursts are denoted as SBl, SB2, and SB3, respectively. 
The star formation rate is assumed to be constantly higher 
for tsbi,s Si t Si Isbi,e iB SBl, tsb2,s St S ^sb2,e Ib SB2, and 
^sb 3 ,s Sts tshz,e iB SB3. The star formation coefficient 
(Csf) is thus described as follows: 


Gsf = 


Cq for quiescent phase 
Csbi for SBl 
ab 2 for SB2 
ab3 forSB3 


( 8 ) 


In the present study, Csf is given for each model (e.g., 
SMC, LMC, and MW). 

We adopt the following formula for the rate of gas ac¬ 
cretion onto the disk of a galaxy: 

A(t) = C'aexp(-t/ta), (9) 

where t^ is a free parameter controlling the timescale of the 
gas accretion and Ca is the normalization factor and deter¬ 
mined such that the total gas mass (Mg_t) accreted onto 
the galaxy can be 1 for a given fa- This Mg,t is determined 
such that Mg^t is the total baryonic mass of a (present) 


galaxy (i.e., gas and stars, including gas , metals, and dust 
ejected from the galaxy through outflow). The parameter 
ta can therefore control the time evolution of the total gas 
mass and the gas mass fraction of a galaxy. We mainly 
show the results of the models with t^ = 5 Gyr, though we 
investigated models with different t^. The main purpose is 
to discuss how dust wind plays key roles in the evolution 
of extinction curves of galaxies in the present study. Such 
roles of dust wind do not depend on the adopted t^. 

Although we show only the results of the models with 
ta = 5 Gyr, the final (i.e., present) gas masses and gas 
mass fractions of the three galaxies are roughly consis¬ 
tent with the observed values. For example, if we adopt 
Mg = 3 X 10 ®Mq (BT12) for the LMG, then the final 
gas mass is 1.1 x IO^Mq and the gas fraction is 0.37. 
These numbers are roughly consistent with the observed 
values (e.g., van den Bergh 2000), which suggests that the 
adopted t^ is reasonable. Furthermore the typical star for¬ 
mation rate at non-starburst phases in the LMC models 
is 0.15 Mq yr“^, which is similar to the observed present 
SFR of the LMC (O.26M0 yr“^; Kennicutt et al. 1995). 
The final gas mass fractions of the SMC and the MW are 
0.53 and 0.05 in the present models, respectively, which 
are also roughly consistent with the observed values (e.g., 
van den Bergh 2000). The evolution of SFRs and gas mass 
fractions are given in Figure 1 and later discussed. 

The assumption of slow gas accretion onto the SMC im¬ 
plies that the SMC has a disk component. Almost all pre¬ 
vious theoretical models of the Magellanic Stream (MS) 
formation needed to assume the gas and stellar disks (e.g., 
Gardiner & Noguchi 1995; Diaz & Bekki 2012) in order 
to reproduce quite well the observed properties of the MS. 
Furthermore, recent observations have confirmed that the 
SMG has a stellar disk (e.g., Dobbie et al. 2014), and it 
is a well known fact that the SMC has a rotating HI disk 
(e.g., Stanmirovic et al. 2004). Therefore, it is reasonable 
for the present study to assume that the SMC has a disk 
onto which gas can be accreted. 


2.4. SN-driven wind 


We assume that some fractions of metals from SNe can 
be removed from a galaxy thorough energetic outflows so 
that they do not contribute to chemical enrichment of the 
galaxy. In the present study only SNe ejecta can be ex¬ 
pelled from a galaxy (“selective wind model”). We do not 
consider that the ejection efficiency depends on (i) chem¬ 
ical elements (e.g., C and O) and (ii) whether the ejecta 
is from SN la or SN II just for simplicity in the present 
study. Some fraction, 1 — /ej, of gaseous ejecta from SNe 
can be mixed with the ISM for chemical enrichment pro¬ 
cesses whereas all AGB ejecta can be mixed with the ISM 
in the selective wind models. 

The SN wind rate {w{t)) is estimated only from the to¬ 
tal mass of gaseous ejecta from SNe (Mej^sn) at each time 
step. The total mass, Mw, of ISM that is removed from a 
galaxy at each time step is as follows: 

Mw = /ejMej,,n. (10) 

Therefore, the SN wind rate at each time step is simply 


w{t) 


dMy„ 

dt 


( 11 ) 


We show the results of the models with /ej = 0.4 in the 
present study, because BT12 have already investigated the 
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SN wind models with /ej = 0.2 and 0.4 and confirmed that 
such models can better reproduce the observed properties 
of chemical abundances of the LMC (e.g., Figure 16 of 
BT12). We also confirm that the roles of dust wind in the 
evolution of extinction curves do not depend strongly on 
/ej, though the models with lower /ej end up with higher 
gaseous metallicities. 

Multiple SN events can strongly influence the chemical 
and dynamical evolution of star-forming regions, and one 
of the most dramatical influences is the formation of su¬ 
per shells and bubbles in the H i disks of the LMC and 
the SMC (e.g., Kim et al. 1999; Stanimirovic et al. 1999). 
Such SN events can cause more efficient dust destruction, 
shattering, and coagulation in shells and bubbles of MCs, 
which implies that the models of dust evolution would 
need to include such dust-related physical processes de¬ 
pending on the time evolution of SN rates. In the present 
study, we do not intend to model such dust-related physi¬ 
cal processes dependent on SN rates, mainly because such 
complicated models require a few additional model param¬ 
eters and thus would possibly prevent us from clarifying 
the key parameters for the evolution of dust extinction 
curves. Thus it would be our future study to include the 
dust-related physical processes dependent on SN rates of 
galaxies in dust evolution models. 

2.5. Delay time distribution of SN la 

We assume that there is a time delay (tia) between the 
star formation and the metal ejection from SNe la. We 
here adopt the following delay time distribution, “DTD”, 
(ff(^ia)) for 0.1 Gyr < tia < 10 Gyr, which is deduced from 
recent observational studies on the SN la rate in extra¬ 
galaxies (e.g., Maoz et al. 2010); 

giaiha) = C'g^Ia > (12) 

where Cg is a normalization constant that is determined 
by the number of SN la per unit mass (which is controlled 
by the IMF and the binary fraction for intermediate-mass 
stars for the adopted power-law slope of —1). We mainly 
investigate models with the above DTD (rather than those 
with typical time delay of SN la being ~ 1 Gyr). The 
fraction of the stars that eventually produce SNe la for 
3-8Mq has not been observationally determined and thus 
is regarded as a free parameter, /b. We investigate models 
with /b = 0.05, because such models can better explain the 
observed chemical properties of the MGs and the Galaxy 
(e.g., Tsujimoto et al. 2010; BT12). 

2.6. Nucleosynthesis and dust yields 

We adopt the nucleosynthesis yields of SNe II and la 
from T95 to deduce j/uy and for a given IMF. For 
AGB yields, we calculate j/agb.i and /i(tagb) for a given 
IMF based on the yield tables from van den Hoek & Groe- 
newegen (1997). Since our recent chemodynamical simu¬ 
lations of disk galaxy formation can reproduce the basic 
dust properties of galaxies (e.g., dust-to-gas-ratios) rea¬ 
sonably well (Bekki 2013), we adopt the dust yields (cor¬ 
responding to /c,i) used in the simulations for the present 
study. The adopted fc,i is originally from Dwek (1998), 
which explained the dust properties of the Galaxy quite 
well. The adopted values are different between different 


elements and between SNe and AGB stars. We investi¬ 
gate the time evolution of dust and metals for the selected 
eight elements, G, O, Mg, Fe, Ga, SI, Ti, and S. 

2.7. IMF 

We adopt a standard power-law IMF; 

4'(toi) = (13) 

where mi is the initial mass of each individual star and the 
slope a is set to be 2.35, which corresponds to the Salpeter 
IMF (Salpeter 1955). The normalization factor Ms^ is a 
function of a, mi (lower mass cut-off), and mu (upper mass 
cut-off). These mj and mu are set to be O.IMq and 5OM0, 
respectively. As discussed in BT12, the IMF slope a can 
influence the chemical evolution history of the LMG. We 
however did not discuss such an influence in the present 
study, firstly because the purpose of this paper is not to 
clarify the roles of IMF in galactic chemical evolution and 
secondly because the models with the standard IMF (with 
SN wind) can explain the observed abundances of the MGs 
(e.g., Tsujimoto & Bekki 2009; BT12). 

The IMF slope can possibly influence the evolution of 
dust and metals of galaxies, because it can influence the to¬ 
tal amount of energy ejected from massive stars and SNe. 
If we adopt a top-heavy IMF (e.g., a = 2.1 instead of 
2.35), then dust and metals can be more efficiently ejected 
from the MGs so that dust-to-gas-ratios can be lower in 
the MGs. The time evolution of extinction curves of the 
MGs might be changed to some extent owing to the more 
efficient removal of dust. Although this is an important is¬ 
sue, we do not investigate how the IMF can influence the 
evolution of dust and extinction curves of the MGs in the 
present study: It is our future study to clarify these IMF 
influences. 

2.8. Two-composition approximation 

Although we investigate the time evolution of small and 
large grains for each element (e.g., G, O, and Mg) sepa¬ 
rately, we divide these grains into the following two cate¬ 
gories: G-grains and Si-grains. Dust grains formed from 
the element G are referred to as G-grains whereas those 
from elements other than G are referred to as Si-grains. 
Asano et al. (2013) also adopted this treatment of Si- 
grains, expecting that solid materials other than carbona¬ 
ceous dust compose silicate materials after dust processing 
in the ISM. Accordingly, there are four different types of 
grains {four-component dust model), for each of which the 
total mass and dust-to-gas-ratio (D) are investigated in 
the present study. For example, the dust-to-gas-ratio for 
small Si-grains is estimated as follows: 

i?s.Si = ^A,i, (14) 

i 

where the summation is done for elements other than G 
(i.e., only for i=0, Fe, Mg, Ga, Si, Ti, and S). The dust- 
to-gas-ratio for small G-grains is simply Ds,c- Thus, the 
total dust-to-gas-ratio [D) is given as: 

^ = ( 15 ) 

i 3 

where summation is done over i=s (small) and 1 (large) 
and j=C, O, Fe, Mg, Ga, Si, Ti, and S. 
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2.9. Dust growth 


We consider (i) that dust grains can grow quite effi¬ 
ciently through the accretion of gas-phase metals onto the 
grains in the cold molecular clouds, (ii) that the timescale 
is shorter for clouds with higher metal-abundance, and (hi) 
that dust accretion timescale (Tacc.s.i and Tacc,i,i) depends 
on the metallicity and the molecular gas density in ISM. 
The accretion timescale for the small grain of ith element 
is therefore as follows: 


_ Tacc,s,2,0 

Tacc,s,2 — 5 

^acc,2 


(16) 


where Tacc.s.i,© is the dust accretion timescale at the solar 
metallicity (0.014) and the total dust abundance (0.0064) 
at the solar neighborhood (in the present MW), and the 
correction factor Caccy is as follows: 


Cacc,* = i^(^z.g,/Hj 


■^Z,g,© 



(17) 


where D, are the dust-to-gas-ratio, the gas-phase 

metallicity, and the cold gas (molecular gas) fraction, re¬ 
spectively, and Zz g 0 and Dq are the gas-phase metallic¬ 
ity and the dust-to-gas-ratio at the solar neighborhood, re¬ 
spectively. Here, Dq is the total dust-to-gas-ratio and thus 
set to be 0.0064, and (3 relates the dust-to-gas-ratio to the 
molecular fraction of ISM. In this formula, we adopt a very 
reasonable assumption that the dust accretion timescale is 
shorter in higher metallicity and higher cold gas fraction. 
The cold gas fraction is assumed to be a function of D as 
follows: 

/h 2 oc (18) 

Uq 


where we estimate /hs by using D and gas mass at each 
time step in a model. We choose this /3 of 0.2 for be¬ 
cause such a choice is consistent with recent observational 
results by Corbelli et al. (2012) who show a correlation 
between the total masses of molecular hydrogen and dust 
in galaxies. 

Although we adopt the exactly the same formula for 
the dust accretion timescale (Tacc,i,i) of large grains, we 
consider that Tacc,i,i,© should be significantly longer than 
Tacc,s,i,© owing to the possible mass growth timescale pro¬ 
portional to dust sizes (H15). We therefore assume that 
Tacc,!.!,© is lOTacc.s.i,© for all models. We show the results 
of the models with Tacc.sy,© = 0.02 Gyr Tacc.iy,© = 0.2 Gyr 
for all elements in the present study. The dust destruction 
timescale (tsn.s.i and tsn,!,*) is set to be 0.2 Gyr for small 
and large grains of all elements, and these values are simi¬ 
lar to those adopted in previous studies (e.g., Dwek 1998; 
H15). 


2.10. Coagulation and shattering 

As discussed by H15, the timescales of coagulation and 
shattering depend on D (dust-to-gas-ratio) in the ISM of 
galaxies. We therefore adopt the following formula: 

Tsh 0 = Tsh,t,©( )~\ (19) 

-^la,© 

where Tshy,© and 7?i,i,© are Tshy and D\^i (dust-to-gas-ratio 
for large grains) at the solar neighborhood, respectively. 
Since C (carbon) becomes carbon dust whereas metals 


other than C become silicate, the shattering timescale for 
large C-grains is as follows: 


^sh,z — ^sh,C,ol y-v ) 5 

M,© 


( 20 ) 


where Di^c is the dust-to-gas-ratio of large G-grains and 
£> 1,0 is the dust-to-gas-ratio for large grains (0.0045) at 
the solar neighborhood. The timescale for elements other 
than G is as follows: 


7sh,i — TshjSi,©! y-, ) ; 


D 


1,© 


( 21 ) 


where A,Si is the dust-to-gas-ratio of large C-grains. We 
choose Tsh,©,© = ©sh,Si =0.1 Gyr, which is similar to those 
adopted in previous studies (e.g., H15). 

Similarly, the coagulation timescale for C is described 
as follows: 


- ('I 

'^co,i — '^cOjC,©! ) 


D 


s,Q 


( 22 ) 


where £>i,c is the dust-to-gas-ratio of small C-grains and 
A,© is the dust-to-gas-ratio for small grains (0.0019) at 
the solar neighborhood. The coagulation timescale for el¬ 
ements other C is as follows: 


_ _ 

©co,i — '^co,Si,©l y-, ) j 


D. 


s,Q 


(23) 


where A.Si is the dust-to-gas-ratio of small Si-grains. The 
values of Tco,©,© and Tco,si,© are set to be 0.02 Gyr, which 
is similar to those adopted by other studies (e.g., H15). 


2.11. Radiation-driven dust wind 

The new model for dust wind rate ITd(t) in the present 
study is different to some extent from those adopted in our 
previous study (BT14). The wind rate is assumed to be 
different not only between C- and Si-grains but also be¬ 
tween small and large grains in the new model. Therefore 
the wind rate Wd(t) is described as follows: 

Wdit) = C^S{t)Md{t), (24) 

where Gw is a parameter that controls the removal effi¬ 
ciency of dust through dust wind, S{t) is the total lumi¬ 
nosity (proportional to the strength of radiation pressure 
of stars), and is the total dust mass. The physical 
meaning of Gw, which is the most important parameter 
in the present study, is clearly explained and discussed 
in Appendix A. About 0.4% of dust can be lost through 
radiation-driven stellar wind within 10 ® yr in a model with 
Gw = 0 . 01 . 

The value of Gw can be different between small and 
large C-grains and Si-grains. For example, the wind rate 
for small C-grains (Gw,s,c; * =C) is given as: 

bFd,s,c(0 = C'w,s,C'S'(t)Md,s.c, (25) 

where Gw,s,c is the value of Gw for small C-grains. As 
demonstrated by recent numerical simulations of dust 
wind in galaxies (Bekki 2015), it depends on a number 
of factors whether dust can be ejected and removed from 
galaxies beyond their halo regions. For example, gas-dust 
interaction through gaseous drag can prevent dust from 
being ejected from gas disks of galaxies. Since this gas- 
dust interaction can not be modeled in the present one- 
zone models, we regard Gw as a free parameter. Thus, we 
investigate the effects of dust wind in the evolution of dust 
by changing the Gw parameters for different grains. 

In the present study, we assume that dust wind can oc¬ 
cur only in strong starburst phases (SBl, SB2, and SB3): 



7 


no dust can be removed from a galaxy in the quiescent 
star-formation. This assumption is quite reasonable and 
realistic, because dust wind is observed in starburst galax¬ 
ies like M82 (e.g., Kaneda et al. 2010). The strength of 
radiation pressure of stars in a galaxy is assumed to be 
proportional to the total luminosity of the galaxy, and the 
total luminosity is calculated as follows: 

Sit) = fM/Lit - T)MUT)dT, (26) 

where /m/l is the mass-to-light ratio {M/L) of a single 
stellar population (SSP) and Mns(T) is the total mass of 
stars formed at t = T. Since the M/L of a SSP is a func¬ 
tion of age and metallicity, we need to use a stellar pop¬ 
ulation synthesis code in order to properly calculate /m/l 
at each time step. Accordingly, we use the code MILES 
(Vazdekis et al. 2010) for the M/L estimation of all mod¬ 
els in the present study. We adopt the SSP table for the 
Salpeter IMF and [Z/H]=—0.4 from the MILES and esti¬ 
mate M/L for stars with different ages in all models. 

Since the parameter Cw is already introduced in equa¬ 
tion (23), the above equation (25) is slightly different from 
that (with a parameter Cr) used in BT14. The radiation- 
driven wind can remove dust from a galaxy, and thus Cw is 
non-zero only for tsbi.s <t< tsbi.e, isb 2 ,s <t< fsb 2 ,e, and 
tsh 3 ,s < t < tsh 3 ,e- C'w = 0 in the quiescent SF phase. We 
have run a large number of models with different values 
of C'w.s.c, C'w.i.Cj C'w,s,Si, and Cw,i,si so that we can find 
models that can explain the observed extinction curves of 
the LMC and the SMC. We mainly show the results of 
three “Hducial” models with a particular combinations of 
these four parameters for the SMC, the LMC, and the 
MW, because they can reproduce the observed extinction 
curves. 


2.12. Derivation of extinction curves 


Based on the total masses (or fractions) of small and 
large C-grains and Si-grains in a model, we can draw 
the extinction curve by using the two-size approximation 
method by H15. Since the details of the method to de¬ 
rive extinction curves from dust masses and compositions 
are already given in H15, we briefly describe it here. We 
adopt the “modified” log-normal model for the grain size 
distribution of small and large grains, and the functional 
form (nij(a)) for the distribution is given as follows: 




= ^^exp{- 


2a2 


}, 


(27) 


where a is the grain size, subscript i describes the small 
(i=s) or large (f=l) grains, subscript j describes carbona¬ 
ceous (/=C) or silicate (j=Si) grains, Cij is the normaliza¬ 
tion constant, and 09 ,^,/ and a are the central grain radius 
and the standard deviation of the log-normal distribution, 
respectively. It should be noted here that the mass distri¬ 
bution (oc a^nij{a)) is log-normal in the adopted model. 

We adopt ao„s,c = ao,.s,Si = 0.005 /im, ao„i,c = ao„i.Si = 
0.1 /im, and a = 0.75 so that the size distribution function 
can roughly cover the small and large grain size ranges. 
The normalizing constant for small and large C-grains and 
Si-grains are determined by the following equation: 


/■°° 4 

limnDij = / -TTa^snij{a)da, (28) 

Jo 3 


where /r = 1.4 is the gas mass per hydrogen nucleus, mn 
is the hydrogen atom mass, and s is the material density 
of a dust grain. Accordingly, we can derive Cij for Dij 
derived from one-zone chemical evolution models by using 
the above equation. 

Based on Uij estimated from Dij at each time step in a 
chemical evolution model, we can derive the dust extinc¬ 
tion curve at the time step. The extinction at wavelength 
A in units of magnitude (A^) normalized to the column 
density of hydrogen nuclei (iVn) is described as follows 
(H15): 

Aa 


= 2.51ogeEE 




where Qextia, X) is the extinction efficiency factor. The 
values of this Qext factor at each wavelength and for each 
grain are evaluated by using the Mie theory (Bohren & 
Huffman 1983) and the same optical constants for silicate 
and carbonaceous dust in Weingartner & Draine (2001). 
We adopt s = 3.5 g cm“^ for small and large Si-grains and 
s = 2.24 g cm“^ for small and large C-grains. As shown in 
H15, this extinction model can reproduce reasonably well 
the observed MW extinction curve. 


2.13. The model parameters for the MCs and the MW 

We first investigated many models with different model 
parameters for dust growth, shattering, and coagulation 
and found that no models can reproduce the observed ex¬ 
tinction curves without the 2175 A bump (SMC) and 
the very weak one (LMC). We then considered the effects 
of dust wind on the evolution of dust extinction curves 
by incorporating the dust ejection/removal process asso¬ 
ciated with dust wind into the one-zone models. Hou & 
Hirashita (2015) have adopted one-zone models withont 
dust wind for deriving dust extinction curves of galaxies 
and found that the models with a wide range of parameters 
for dust-related physical processes (e.g., coagulation) can 
not reproduce the observed extinction curves without the 
2175 A bump. Therefore, it seems that standard chemical 
evolution models without dust wind can not reproduce the 
extinction curves of the MCs. 

After extensive investigation of models with dust-related 
physical processes (e.g., selective destruction of C-grains 
by SNe and dust wind), we could finally find that models 
with dust wind can reproduce better the observed extinc¬ 
tion curves of the MCs. Therefore, we focus on the results 
of the models with dust wind in the present study. Al¬ 
though it is possible that physical processes other than 
dust wind can be responsible for the origin of the extinc¬ 
tion curves of the MCs, we focus exclusively on the effects 
of dust wind on the evolution of extinction curves of the 
MCs. 

We consider that the efficient ejection/removal of dust 
grains from the SMC and the LMC can be the key phys¬ 
ical process for the origin of their dust extinction curves. 
We therefore (i) change the parameters of dust wind and 
(ii) adopt reasonable values for other model parameters 
(e.g., for dust growth) so that we can more clearly elu¬ 
cidate the key factor for the presence or absence of the 
2175 A bump. The model parameters for star forma¬ 
tion (e.g., burst epoch) and dust-related physical processes 
(e.g., dust accretion timescale) are summarized in Table 1 



and Table 2, respectively, for the MCs and the MW. BT12 
have already discussed how the chemical evolution (thus 
dust evolution) histories of the MCs depend on star for¬ 
mation parameters, and H15 has discussed how the dust 
evolution of a galaxy depends on the model parameters 
for dust growth, destruction, shattering and coagulation. 
Furthermore, we have confirmed that only some models 
with dust wind can explain the observed extinction curves 
of the MCs. Therefore, we focus exclusively on how dust 
wind plays a key role in controlling the extinction curves 
of the MCs. 

In determining the three burst epochs, we consider both 
recent observational and theoretical results on the star 
formation and chemical evolution histories of the MCs. 
Harris & Zaritsky (2004) investigated the ages of stars 
in the SMC based on UBVI photometry from their Mag¬ 
ellanic Clouds Photometric Survey and found signihcant 
rises in the mean star formation rate around 0.06, 0.4, 
and 2.5 Gyr ago. Using gas dynamical simulations of the 
SMC interacting with the LMC and the MW, Yoshizawa & 
Noguchi (2003) predicted that the SMC experienced star- 
burst events about 0.2 Gyr and 1.5 Gyr ago. Rubele et al. 
(2015) have recently detected the starburst populations 
with ages of 1.5 and 5 Gyr in the deep images of the SMC 
from the VISTA survey of the Magellanic Clouds in the 
YJKg) bands. Guided by these observations, we assume 
that the SMG experienced three enhanced star formation 
episodes at Gb.i.s = 8.5, 11.5, and 12.85 Gyr (assuming 
the current age of 13 Gyr for the LMG). 

The LMG also appears to have experienced at least 
three major epochs of significantly enhanced star forma¬ 
tion. Rubele et al. (2012) investigated the AMR of the 
LMG stellar populations and found that the SFR peak 
around 2 Gyr ago, though the AMR appears to be different 
between different local regions. Based on the results of the 
gas dynamical simulations of the LMC over the past 6 Gyr, 
Bekki et al. (2004) suggested that a starburst, which can 
form globular clusters, could have occurred in the LMC 
about 3 Gyr ago owing to the LMC-SMC-Galaxy tidal in¬ 
teraction. Harris & Zaritsky (2009) found very strong en¬ 
hancement of mean star formation from 10® yr ago to the 
present in the blue arm. Constellation HI, and 30 Doradus 
region of the LMC. Thus, we assume that the LMC experi¬ 
enced three enhanced star formation at tsb,i,s = 10.0, 11.5, 
and 12.85 Gyr (corresponding to 3, 1.5, and 0.15 Gyr ago, 
respectively). 

Table 3 summarizes the model parameters of dust wind 
in the three galaxies. Although we investigated much more 
models than those listed in this table, we show only these 
representative ones, because they grasp some essential in¬ 
gredients of the dust wind effects on the extinction curves 
of galaxies. We mainly describe the time evolution of the 
extinction curve in the fiducial model for the SMC (SI). 
This fiducial model is chosen, because it can best repro¬ 
duce the observed extinction curve of the SMC. As a com¬ 
parative experiment, the model without dust wind (yet 
with starburst; S2) is investigated for the SMC. These 
models with and without dust wind enable us to demon¬ 
strate more clearly how the dust wind can change the 
shapes of the extinction curves of the SMC. The models 
with different C'w,s,c (S3, S4, and S5) are investigated so 
that the importance of ejection/removal of small C-grains 


in the evolution of the extinction curve of the SMC can 
be more clearly described. The models without starbursts 
(S6—S9) enables us to understand the importance of star- 
burst dust wind in the evolution of the extinction curve. 

We show the results of the LMC models with and with¬ 
out dust wind also in order to clarify the physical rea¬ 
sons for the observed weak 2175 A extinction bump in 
the LMC. We show the results of the MW model (Ml) 
without starburst (thus without dust wind) only, because 
there is little observational evidence for recent strong star- 
burst in the MW. This MW model and those for the MCs 
(SI and Ml) can combine to demonstrate that the present 
models can explain not only the extinction curve of the 
MW but also those of the MCs. Figure 1 shows the time 
evolution of Mg (normalized), /g (gas mass fraction), and 
SFR (normalized) for the fiducial SMC, LMC, and MW 
models. 

3. RESULTS 
3.1. SMC 

Figures 2 and 3 show the time evolution of the extinc¬ 
tion curve and the relative fractions of dust grains, respec¬ 
tively, in the standard SMC model (SI) in which the SMC 
can experience three strong starburst (SB) events. Before 
the first SB event (SBl), the star formation rate of the 
SMC can be kept low so that chemical evolution can pro¬ 
ceed very slowly (T < 8.5 Gyr). In this low-metallicity 
phase (pre-SBl phase), stellar dust production dominates 
the dust content (i.e., yet little dust growth) since the dust 
abundance is too low for other dust processing mechanisms 
(in particular, shattering) to work efficiently. As a result 
of this, the ISM dust of the SMG is dominated by large 
grains. The SMC thus can show a rather flat extinction 
curve without the 2175 A bump in this pre-SBl phase. 

As the gas-phase metallicity becomes higher through 
more rapid chemical evolution during SBl, the mass frac¬ 
tion of small grains becomes larger owing to the interplay 
between shattering and accretion (T = 9.8 Gyr). The ex¬ 
tinction curve can therefore show a slightly steeper rise 
in the FUV regime and the 2175 A bump can be appre¬ 
ciably seen at this post-SB 1 phase. The extinction curve 
becomes progressively steeper and the 2175 A bump be¬ 
come more prominent after the SBl (T = 11.5 Gyr). The 
fraction of cold gas can become higher at this epoch, be¬ 
cause D, which can determine the gas fraction is higher. 
The fraction of small Si-grains can be as large as that of 
large Si-grains at this pre-SB2 phase, because the higher 
metallicity and the higher cold gas fraction combine to al¬ 
low the more rapid growth of the small grain abundance. 
These behaviors in the evolution of Si grains are appar¬ 
ently similar to those reported in Asano et al. (2014). 

Although the fraction of small Si-grains can be tem¬ 
porarily lower owing to the efficient destruction of dust 
by SNe during the second starburst phase (SB2), it can 
start to rise again soon after SB2 (T = 12.5 Gyr). The 
total mass of small Si-grains becomes 1.6 times larger than 
that of large Si-grains at T = 12.5 Gyr so that the mass 
fraction of small Si-grains can be the largest in the SMG 
history. Gonsequently, the extinction curve of the SMG in 
this model at T = 12.5 Gyr shows a very steep rise that is 
even steeper than the observed SMG extinction curve and 
thus inconsistent with the observed curve. Although small 
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C-grains can also suffer a significant loss of their mass ow¬ 
ing to the grain destruction by SNe during SB2, their mass 
can steadily and rapidly increase through the accretion of 
gas-phase metals after SB2. This increase of small C-grain 
abundance is responsible for the strong 2175 A bump in 
the extinction curve at T = 12.5 Gyr. 

The fraction of small Si-grains can become smaller ow¬ 
ing to their destruction by SNe and ejection from the SMC 
through radiation-driven dust wind in the third starburst 
phase (SB3) so that the extinction curve at T = 12.9 Gyr 
can be slightly less steeper than that at T = 12.5 Gyr. The 
dust wind effects are more significant in this SB3, mainly 
because the total luminosity of stars in the SMC takes 
its maximum. Although a significant fraction of small C- 
grains can be removed from the SMC through dust wind, 
the extinction curve still clearly shows the 2175 A bump 
in the early phase of this SB3. The 2175 A bump can com¬ 
pletely disappear only after most of the small C-grains can 
be removed from the SMC through dust wind at T = 12.9 
Gyr. The final extinction curve at T = 13 Gyr is very 
similar to the observed one both in the lack of the 2175 
A bump and in the steep linear rise at FUV wavelength. 
Thus the present standard SMC model with dust wind can 
reproduce quite well the observed extinction curve. 

As clearly shown in Figure 3, the time evolution of C- 
grains is different between the SMC models with and with¬ 
out dust wind (SI and S2, respectively) after SB2: the 
fractions of small and large C-grains and small Si-grains 
(relative to large Si-grains) are not so different between 
the two models during and after SBl and SB2. This sug¬ 
gests that the SMC can not show the extinction curve 
without the 2175 A bump until quite recently. This also 
suggests (i) that star-bursting dwarf galaxies do not nec¬ 
essarily show extinction curves without the 2175 A bump 
and (ii) that dust wind needs to be strong enough to re¬ 
move most (90 percent) of small C-grains in star-bursting 
dwarfs with no 2175 A bump in their extinction curves. 
It is thus highly likely that starburst dwarfs have different 
strength of the 2175 A bump. 

Figure 4 shows that the final extinction curve in the 
SMC model with three SB events yet no dust wind has 
the conspicuous 2175 A bump. The SMC models with 
lower dust wind efficiency (C'w,s,c=0.01 and 0.03 for S4 
and S5, respectively) also show the 2175 A bump and 
the FUV rise that is less consistent with the observed ex¬ 
tinction curve of the SMC. The SMC model with higher 
dust wind efficiency (C'w,s,c=0-1, S6) can reproduce the 
observed extinction curve of the SMC as well as the stan¬ 
dard SMC model. These results confirm that the removal 
of small C-grains from the SMC through dust wind is a vi¬ 
able mechanism for the observed lack of the 2175 A bump 
in the extinction curve of the SMC. 

The model with no starburst thus no dust wind (S7) 
shows the conspicuous 2175 A bump and a steeper FUV 
rise, which clearly demonstrates that dust wind can signif¬ 
icantly influence the shape of the extinction curve of the 
SMC. As shown in the models S8 and S9, the presence or 
absence of the 2175 A bump dose not depend on whether 
SBl or SB2 occurs in the SMC. However, the 2175 A 
bump of the extinction curve is too strong in the model 
SIO without SB3, which suggests that the strength of the 
2175 A bump can depend on whether the SMC experi¬ 


ences SB3 with dust wind. These results thus confirm that 
the last SB with strong dust wind about 0.2 Gyr ago is a 
key physical process for the observed characteristic extinc¬ 
tion curve of the SMC. The results of the SMC models with 
different strengths of starbursts are discussed in Appendix 
B. The metallicity evolution of the SMC (LMC and MW) 
in the present new models with dust evolution is discussed 
briefly in Appendix C. 

3.2. LMC 

Figure 5 presents the time evolution of the extinction 
curve of the LMC in the model LI with three SB events. 
Gas can be more rapidly converted into new stars in this 
LI than in the SMC models (owing to the adopted higher 
SF efficiency) so that chemical abundances can be higher 
even well before the first SB event (T = 8.4 Gyr). As 
shown in Figure 6, small C- and Si-grains can start to 
grow earlier in this LI, and their mass fraction can take 
their peak values before SBl. Consequently, the extinc¬ 
tion curve can have a rather conspicuous 2175 A bump 
and a very steep FUV rise at T = 8.4 and 9.8 Gyr. Since 
SNe can efficiently destroy dust grains during SBl owing 
to high SF rates (i.e. high SN rates), the 2175 A bump 
becomes weakened after SBl (T = 11.5 Gyr). However, 
such destruction of small grains can not severely weaken 
the 2175 A bump after SBl and even after SB2 (T = 12.5 
Gyr). 

The 2175 A bump of the LMC becomes less conspicu¬ 
ous during SB3 owing to the combination effect of dust de¬ 
struction and radiation-driven dust wind (T = 12.9 Gyr). 
The shape of the extinction curve of the LMC can not be¬ 
come quite similar to the observed one until the end of SB3 
{T = 13 Gyr). The rather weak 2175 A bump at T = 13 
Gyr in this model is due to the more preferential removal 
of small C-grains from the LMC through radiation-driven 
dust wind during the longer starburst (SB3). The ma¬ 
jor difference in the dust abundances between this LMC 
model LI and the SMC model SI is that the final frac¬ 
tions of small and large C-grains are both slightly larger 
in LI than in SI (See Figures 4 and 6). These differences 
can cause a difference in the extinction curves between the 
LMC and the SMC in the present study. 

Figure 7 compares between the extinction curve of LI 
with dust wind and that of L2 without dust wind. Clearly, 
LI with dust wind can better reproduce the observed ex¬ 
tinction curve of the LMC, which confirms that the effi¬ 
cient removal of C-grains through radiation driven dust 
wind is essential for explaining the origin of the extinc¬ 
tion curve of the LMC. Previous numerical simulations of 
the MCs (e.g., Yoshizawa & Noguchi 2003; Bekki & Chiba 
2005) suggested that the observed enhanced star forma¬ 
tion rates in the MCs (e.g., Harris & Zaritsky 2006, 2009) 
can be due largely to the LMC-SMC-Galaxy tidal inter¬ 
action in the past 3 Gyr. Thus the results in Figures 2-7 
suggest that the origin of the observed extinction curves of 
the LMC and the SMC can be closely associated with the 
past active star formation triggered by tidal interaction 
among the three galaxies. 

3.3. MW 

It is instructive for the present study to show whether 
the new models for dust evolution of galaxies can repro- 
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duce the observed extinction curve of the MW just by 
changing the model parameters. Figure 8 describes the 
time evolution of the MW model Ml without starburst 
events (thus without dust wind). In the very early chem¬ 
ical evolution history of the model Ml (T = 1 Gyr), the 
extinction curve is rather flat without the clear 2175 A 
bump owing to the very small fraction of small C-grains. 
The extinction curve of the model Ml evolves to have the 
2175 A bump and a steeper rise at 1/A > 6 /im“^ only 
2 Gyr after the commencement of its star formation (i.e., 
disk formation). The overall feature of the extinction curve 
becomes quite similar to the observed one for the MW at 
T = 5 Gyr, which implies that MW-like galaxies at higher 
redshifts might already have extinction curves similar to 
that of the present MW. 

The 2175 A bump is a bit too conspicuous at T = 5 Gyr 
in comparison with the observed one, and the bump does 
not evolve significantly with time after T = 5 Gyr. This 
is mainly because the relative mass fractions of small Si- 
grains and small and large G-grains (with respect to large 
Si-grains) do not change significantly after T = 5 Gyr. 
The overall shape of the extinction curve at T = 13 Gyr 
becomes closer to that of the observed one. However, the 
final shape of the 2175 A bump is still slightly different 
from the observe one (e.g., more remarkable bump), which 
implies that the present models do not reproduce the dust 
abundances of the MW completely well. It is our future 
study to reproduce better the observed extinction curve 
of the MW by using a more sophisticated dust evolution 
model. 

These results imply that luminous disk galaxies like the 
MW are likely to show clear 2175 A bumps in their ex¬ 
tinction curves, because a significant amount of dust is 
unlikely to be removed from the galaxies owing to their 
deep gravitational potential. Although a small amount of 
dust and metals could be removed/ejected from the galax¬ 
ies through SN and dust wind, such a less significant re¬ 
moval/ejection process can not influence the shapes of the 
extinction curves to a large extent. The present study 
therefore predicts that the shapes of extinction curves de¬ 
pend on the masses of galaxies in such a way that more 
massive galaxies are likely to have the 2175 A bump. It 
is an interesting observational question how the strength 
of the 2175 A bump in a massive galaxy depends on its 
physical parameter (e.g., total mass and luminosity). It is 
worth noting here that Asano et al. (2014) and Nozawa 
et al. (2015) attempted to reproduce the 2175 A bump of 
the MW by freely changing the mass fractions of various 
ISM phases. 

Finally, Figure 9 summarizes the time evolution of mass 
fractions of the key three dust grains with respect to large 
Si-grains for the MGs and the MW. The differences in 
the evolution between the three galaxies can be clearly 
seen, which reflects on the time evolution of the extinction 
curves of the three. Star formation histories with or with¬ 
out starbursts and capabilities to retain dust within galax¬ 
ies are key determinant for the time evolution of the ex¬ 
tinction curves of galaxies. Since the present study has in¬ 
vestigated only three different types of galaxies (i.e., SMC, 
LMC, and MW), we will need to investigate models with 
different star formation and dust evolution histories to un¬ 
derstand fully the origin of various extinction curves ob¬ 


served in galaxies with different masses and Hubble types. 

4. DISCUSSION 

4.1. The origin of the 2175 A bump 

We have, for the first time, demonstrated that the ob¬ 
served extinction curve without the 2175 A bump in the 
SMC can be reproduced, if the small C-grains can be more 
preferentially lost through radiation-driven wind during 
the last starburst around 0.2 Gyr ago. In this “lost small 
C-grain” scenario, the LMC can have a weak 2175 A 
bump in its dust extinction curve, because it can retain a 
larger amount of small C-grains owing to its deeper gravi¬ 
tational potential in comparison with the SMC. This “lost 
small C-grain” scenario predicts that more massive galax¬ 
ies (like the MW) are likely to show more distinctive 2175 
A bumps in their dust extinction curves. It should be 
stressed, however, that the present low-mass galaxies can 
show such 2175 A bumps if they have not experienced 
recent strong starbursts and the resultant loss of small C- 
grains. 

The key physical process required for the lack of the 
2175 A bump in a galaxy is the more efficient removal 
of C-grains than Si-grains from the galaxy. The removal 
(or stripping) process of dust grains in galaxies should be 
rather complicated, because gravitational dynamics, gas- 
dust drag, and radiation pressure of stars can be all in¬ 
volved in the removal process. Therefore, it would be a 
formidable task to estimate accurately the dust removal 
efficiency in a galaxy for the physical parameters of the 
galaxy. Thus we investigate whether and in what physi¬ 
cal conditions the above two requirements can be met in 
a galaxy by using a simplified yet reasonable model. 

The first question is whether these C-grains can be 
more efficiently removed from a galaxy through radiation- 
driven wind than Si-grains. The key parameter for this 
removal process in a galaxy is the frequency-averaged ra¬ 
diation pressure coefficient (Qpr), as demonstrated in pre¬ 
vious theoretical models and recent numerical simulations 
of dust wind in galaxies (e.g., Ferrara et al. 1991; Blanch & 
Ferrara 2005; Bekki 2015). We have therefore investigated 
Qpj. for C and Si grains by using the publicly available data 
for optical properties of dust (Draine and Lee 1984; Laor 
& Draine 1993) and the stellar population synthesis code 
(“MILES”) developed by Vazdekis et al. (2010). The ra¬ 
diation pressure coefficient (Qpr) for a grain with a radius 
is given as follows: 

Qpr = Qabs + (1 ~ gi^AQsca, (30) 

where Qabs and Qsca are the absorption and scattering co- 
efhcient, respectively, and g{0) is the scattering asymmetry 
parameter. 

Using Qpr and the flux of a SED (/(A)), we can estimate 
Qp^ as follows: 

/ •^max 

QpJ(A)dA, (31) 

min 

where Amin and Amax are set to be 3340.5 A and 7409.6 
A, respectively, in MILES. The above wavelength range 
is publicly available one in MILES and /(A) is the nor¬ 
malized flux (i.e., /(A) = 1). Figure 10 shows 

the wavelength-dependence of radiation pressure coeffi¬ 
cient (Qpr) for silicate and graphite with the radii of 0.01 
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fxm and 0.1 /im and the SED of a galaxy for a metallic- 
ity ([Z/H]) of —0.7 dex and an age of 63 Myr (which is 
the youngest stellar population that can be adopted in 
MILES). We use this SED for a single stellar population 
(SSP) to mimic the SMC with young starburst compo¬ 
nents. We discuss just this model with two representative 
dust radii as an example, though, ideally speaking, we 
should show the results for all dust radii. 

Clearly, (5*^. is higher in graphite than in silicate for the 
adopted SED and dust radii (0.01 /xm and 0.1 /xm), which 
means that graphite is more likely to be influenced by ra¬ 
diation pressure of stars and thus removed from a galaxy 
if the gravitational potential is sufficiently shallow. Fur¬ 
thermore, the difference in (5*j. is significantly larger in the 
small grains (a = 0.01 /xm) than in the large ones (a = 0.1 
^m), which implies that small C-grains can be much more 
efficiently removed than small Si-grains. Although the 
real SED of a SMC-type galaxy might be different from 
the adopted SED in this investigation, it is unlikely that 
this result changes qualitatively if a more realistic SED is 
adopted, given Qp^ of graphite always higher than those 
of silicate for the relevant wavelength range. These results 
thus suggest that graphite is more likely to be lost from 
a SMC-type galaxy than silicate, though a full investiga¬ 
tion on this dust removal process for different dust grains 
and more realistic galaxy spectra needs to be done in our 
future studies. 

If a rather small mass fraction of small C-grains is re¬ 
sponsible for the observed lack of the 2175 A bump in the 
SMC, then is there any observational evidence for that ? 
Using ISOPHOTO and HiRes IRAS data for dust emis¬ 
sion and ATCA data for HI column gas density. Hot et 
al. (2004) investigated the gas-to-dust-ratios of polycyclic 
aromatic hydrocarbon (PAH) molecules or small grains 
(VSGs) for different regions of the SMC. They found that 
the derived difference in the gas-to-dust-ratio of PAHs and 
very small grains between the Galaxy and the SMC is three 
times larger than the metallicity difference between the 
two galaxies. This observational result implies that the 
dust depletion levels for VSG and PAHs are low in the 
SMC for some physical reasons (e.g., destruction by SNe; 
Hot et al. 2004). Since this observational study can not 
distinguish small C-grains from small Si-grains, it is not 
clear whether the observed low depletion level of VSGs 
implies that the mass fraction of small C-grains is signifi¬ 
cantly low owing to efficient removal or destruction of such 
grains. It appears to be difficult for photometric study of 
the SMC alone to estimate the contributions of small C- 
and Si-grains to the observed dusty SED separately. 

The present study predicts that the star-forming regions 
of the SMC for which there is no 2175 A bump can be 
dominated by Si-grains (i.e., small fractions of C-grains). 
However, Welty et al. (2001) found that the Si appears 
to be undepleted in a local region of the SMC toward the 
SMC star Sk 155, which means that the local region should 
contain little Si dust. If this is not just for the Sk 155 re¬ 
gion but for other star-forming regions generally in the 
SMC, then the present scenario and other dust models for 
the SMC like Pei (1992) appears to be ruled out. Thus 
it would be observationally important to make a robust 
estimation of the dust depletion level of Si in many star¬ 
forming regions with and without the 2175 A bump for 


the SMC. 

4.2. Alternative scenarios 

A number of possible carriers of the 2175 A bump have 
been extensively discussed by many authors so far (e.g., 
Whittet 2003). Mishra & Li (2015) have recently found 
a strong correlation between the strength of the 2175 A 
bump and the carbon depletion level ([C/H]dust) in the 
MW and thus suggested that the carrier of the observed 
2175 A bump is graphite or PAH molecules. Their re¬ 
sults are consistent with the present “lost small C-grains” 
scenario, though they did not discuss small and large C- 
grains separately in the context of the origin of the 2175 
A bump. The results imply that the time evolution of 
C-grains in a galaxy can be reflected on the strength of 
the 2175 A bump in the extinction curve of the galaxy. 

The severely dehcient small C-grain of the SMC ISM 
discussed in the present study would be just one of the 
promising ways that explain the observed lack of the 2175 
A bump in the SMC. Nozawa et al. (2015) have recently 
reproduced the observed dust extinction curve without the 
2175 A bump in SDSS J1048-I-4637 by adopting the opti¬ 
cal constant of amorphous carbon for carbonaceous grains 
(See Figure 4 in their paper). Although their dust extinc¬ 
tion curve (Nozawa et al. 2015) is quite different from 
that of the SMC, it would be possible that amorphous 
carbon dust might be responsible for the lack of the 2175 
A bump in the SMC extinction curve. Recent observa¬ 
tional study of dust in the LMC by Herschel (Meixner et 
al. 2010) has shown that their dust model with amor¬ 
phous carbon can better explain their observations for a 
more reasonable dust-to-gas-ratio. It appears to be cur¬ 
rently difficult to model the evolution of amorphous carbon 
by considering the possible physical effects (e.g., radiation 
Helds and gaseous shock) related to the transformation of 
carbon grain properties. 

As demonstrated in the present study, a way to repro¬ 
duce the observed lack of the 2175 A bump in the SMC 
extinction curve is to reduce the relative contribution of 
small C-grains. Therefore, it does not make much differ¬ 
ence whether small C-grains can be lost either through 
energetic dust wind or through some destruction process 
of dust grains. If small C-grains can be more preferen¬ 
tially destroyed or transformed by SNe or gaseous shock 
or strong radiation field in the ISM of the SMC, the 
2175 A bump might not be observed in the SMC. It 
would be possible that the last strong LMC-SMC colli¬ 
sion/interaction about 0.2 Gyr ago (e.g., Diaz & Bekki for 
the latest model), which can strongly compress the ISM of 
the SMC (e.g., Bekki & Chiba 2007), could preferentially 
destroy the small C-grains. It is our future study to inves¬ 
tigate whether small C-grains can be more preferentially 
destroyed by these physical processes in the ISM by using 
pc-scale numerical simulations of the ISM. 

4.3. Incorporation of four-component grain models into 
numerical simulations of galaxy formation and 
evolution 

Interstellar dust can have significant influences on the 
physical processes of galaxy formation and evolution, 
such as the formation of molecular hydrogen in the ISM 
of galaxies (e.g., Gould & Salpeter 1963) and photo- 
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electric heating of the ISM (e.g., Watson 1972). Al¬ 
though these influences have been recently included in 
Nbody-|-hydrodynamical simulations of galaxy formation 
and evolution in a self-consistent manner (Bekki 2013, 
2014, and 2015), the evolution of dust sizes, which can 
control dust-related physical processes of ISM (e.g., Yama- 
sawa et al. 2011), has not been incorporated properly into 
numerical simulations of galaxy formation and evolution. 
The present study has demonstrated that the evolution of 
dust sizes and compositions is a key factor for the evolu¬ 
tion of extinction curves in galaxies. We thus suggest that 
if future numerical simulations of galaxy formation and 
evolution can include the time evolution of dust sizes and 
compositions, such numerical studies can have the follow¬ 
ing improvements in their predictive power. 

First, such simulations with the evolution of dust sizes 
and compositions can predict the SEDs of galaxies in a 
more consistent manner. The SEDs of dusty galaxies are 
strongly influenced by the adopted extinction curves, and 
the shapes of the extinction curves of galaxies depend on 
their star formation histories that control their dust evolu¬ 
tion histories. Therefore, previous SED models of galaxies 
with a fixed extinction curve (e.g., an adoption of the MW 
dust extinction curve) might be much less self-consistent 
and realistic. Construction of SEDs from age and metallic- 
ity distributions of stars and spatial distributions of stars 
in numerical simulations has been done by many authors 
(e.g., Bekki et al. 1999; Bekki Shioya 2001; Johnson 2006; 
Yajima et al. 2014), though they had to assume a constant 
dust-to-metal-ratio and a fixed dust extinction curve. It 
is doubtlessly worthwhile for these simulations to include 
the evolution of dust sizes and compositions for a more 
consistent prediction of galactic SEDs. 

Second, future simulations with dust size evolution can 
provide more accurate predictions of molecular hydrogen 
abundances in the ISM of galaxies. As shown in sophisti¬ 
cated one-zone models of galaxy formation with dust size 
evolution by Yamasawa et al. (2012), the formation effi¬ 
ciency of molecular hydrogen on dust grains in a forming 
galaxy depends on the evolution of the dust size distribu¬ 
tion of the galaxy: since the molecular hydrogen forma¬ 
tion rate on the surface of a dust grain depends on the 
surface area, the net formation rate of molecular hydrogen 
on dust grains in a galaxy relates to the dust extinction 
curve. In order for numerical simulations to predict the 
dust size distribution of a galaxy, they need to follow the 
time evolution of dust grains with different sizes for each 
gas particle. This is a very formidable task, because such 
dust grain evolution is computationally heavy even in one- 
zone chemical evolution models (e.g., Asano et al. 2014). 
A wise way to implement the dust size evolution in numer¬ 
ical simulations of galaxy formation would be to adopt a 
two-size approximation as done in this study. 

Third, the effects of photo-electric heating of dust on 
the thermal evolution of ISM in galaxies can be more 
self-consistently investigated in numerical simulations with 
dust size evolution. The photo-electric heating rate de¬ 
pends on the dust size distribution (not just dust abun¬ 
dances), and the photo-electric heating has been demon¬ 
strated to influence galaxy-wide star formation through its 
influences on thermal history of galactic ISM (e.g., Bekki 
2015). Therefore photo-electric heating would be one of 


the key ingredients that future numerical simulations, in 
particular, those for the evolution of galaxy-wide star for¬ 
mation, should include. Thus, we suggest that incorpo¬ 
ration of dust size evolution into numerical simulations of 
galaxy formation and evolution is essential for understand¬ 
ing various aspects of galaxy formation and evolution. 

5. CONCLUSIONS 

We have constructed a new dust evolution model for the 
MCs by incorporating the two-size (i.e., small and large) 
approximation model for dust size distributions (HIS) into 
one-zone chemical evolution models with dust winds. The 
time evolution of small and large dust grains has been 
self-consistently investigated for each element (e.g., C, O, 
and Si) by considering shattering and coagulation of these 
grains. The extinction curve of a galaxy at each time step 
in a one-zone chemical evolution model has been derived 
from the relative abundances of small and large carbona¬ 
ceous grains (‘C-grains’) and small and large silicate grains 
(‘Si-grains’) by using the optical constants for these grains: 
the galaxy has four-component dust in the ISM. We have 
searched for the model parameters with which the ob¬ 
served extinction curves of the MCs can be successfully 
reproduced. However, we mainly showed the results of 
the models with different wind parameters to demonstrate 
more clearly the key roles of dust wind in the evolution of 
dust extinction curves. The main results are as follows. 

(1) The SMC extinction curve with almost no 2175 A 
bump and a rather steep EUV rise can be reproduced well 
by the present models with SFHs consistent with the ob¬ 
served ones (‘SMC model’), if the small C-grains can be 
more efhciently ejected from the SMC through radiation- 
driven dust wind during the latest starburst about 0.2 Gyr 
ago. Without this selective loss of small C-grains, the mod¬ 
els show conspicuous 2175 A feature regardless of model 
parameters for SFHs and dust physics. These results sug¬ 
gest that the significantly enhanced star formation possi¬ 
bly triggered by the last LMC-SMC tidal interaction can 
be responsible for the origin of the 2175 A bump in the 
SMC. 

(2) The EMC extinction curve with weak 2175 A bump 
and a steeper FUV rise can be also reproduced by the 
present EMC models, if the small C-grains can be ejected 
from the EMC more efficiently than other dust compo¬ 
nents. The dust removal efficiency of small C-grains re¬ 
quired to explain the EMC extinction curve should be by 
a factor of ~ 3 smaller than that required to explain the 
SMC extinction curve. Such a less efficient ejection of 
small C-grains is possible for the EMC owing to the deeper 
gravitational potential well. 

(3) The present model can reproduce the extinction 
curve of the MW with the 2175 A bump if there is no re¬ 
cent starburst event (as observed) that can trigger strong 
dust wind in the model. This is simply because small 
C-grains can not be preferentially lost in the MW model 
with the quiescent star formation history. We have sug¬ 
gested (i) that dust grains are less likely to be ejected from 
the MW owing to its deep gravitational potential well (in 
comparison with the MCs) and (ii) that this more efficient 
trapping of dust grains in the MW can be responsible for 
the evolution of the extinction curve. 

(4) The selective loss of small C-grains from the ISM of 
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the MCs is demonstrated to be physically reasonable and 
realistic, because the frequency-averaged radiation pres¬ 
sure efficiency (Q*!-) is significantly larger in small C-grains 
than in small Si-grains for reasonable spectral energy dis¬ 
tributions of the MCs. Since there are other alternative 
mechanisms for the origin of the extinction curves of the 
MCs with no/little 2175 A bump (e.g., selective destruc¬ 
tion of small C-grains and amorphous carbon grains in the 
MCs), these alternative scenarios will need to be investi¬ 
gated in our future papers. 

(5) The present study suggests that the dust extinction 
curves could be different between galaxies with different 
masses (i.e., difference in the depth of gravitational po¬ 
tential) and different star formation histories (e.g., with 
or without starbursts), because the evolution of dust com¬ 
position and size distributions can be influenced by these 
physical properties of galaxies. We will need a full numer¬ 
ical simulation including most of the dust-related physical 
processes in order to understand the dependences of the 
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dust extinction curves on the physical properties of galax¬ 
ies. 

The strength of the 2175 A bump, their correlations 
with galaxy properties, and the influence of the bump on 
galactic SEDs have been investigated for galaxies with dif¬ 
ferent masses and types at low and high redshifts (e.g., 
Calzetti et al. 1994; Noll et al. 2009; Inoue et al. 2006; 
Conroy 2009; Buat et al. 2011; Wild et al. 2011; Kriek 
& Conroy 2013; Mao et al. 2014; Scoville et al. 2015). 
Although the present study did not discuss these obser¬ 
vations, the new results based on the present dust wind 
model can provide some clues to the origin of these ob¬ 
served properties. We will discuss these issues in our forth¬ 
coming papers. 
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Table 1 

The model parameters for SF histories of the three galaxies. 


Galaxy name 

Gq " 

b 

^sbl 

Gsb2 

d 

^sb3 

^sbl,s 

^sbl,e 

^sb2,s ^ 

^sb3,e 

^sb3,s 

GbS.e 

SMC 

0.005 

0.05 

0.05 

0.1 

8.5 

8.6 

11.5 

11.6 

12.85 

12.99 

LMC 

0.01 

0.05 

0.05 

0.1 

10.0 

10.1 

11.5 

11.6 

12.85 

12.99 

MW 

0.06 

- 

- 

- 

- 

- 

- 

- 

- 

- 


^'The star formation coefficient (Csf) during the quiescent phase (i.e., without starburst). 

*^(781 during the first starburst phase (SBl). The mark in the MW model means Csi is always Cq 
owing to no starburst events. 

^Csf during the second SB phase (SB2). 

'^C'sf during the third starburst phase (SB3). 

®The time (Gyr) when the first SB starts. The mark in the MW model means no SB events in its 
star formation history. 

^The time (Gyr) when the first SB ends. 

®The time (Gyr) when the second SB starts. 

*'The time (Gyr) when the second SB ends. 

‘The time (Gyr) when the third SB starts. 

^The time (Gyr) when the third SB ends. 
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Table 2 

The model parameters for the timescales for the investigated physical processes of dust. 


Physical process 

Parameter ^ 

small Si ’’ 

large Si 

small C 

large C 

Accretion 

"^acc,© 

0.02 

0.2 

0.02 

0.2 

Destruction 

'TSN,© 

0.2 

0.2 

0.2 

0.2 

Shattering 

Tsh,© 

- 

0.1 

- 

0.1 

Coagulation 

Tco,© 

0.2 

- 

0.2 

- 


simple version of the parameter description is given (e.g., race instead of 
Ta.cc,i,j)- The value for each parameter (e.g., Tacc,©) is for the metal and dust 
abundances of the solar neighborhood (in the present Galaxy). 

’’Dust is assumed to consist of four different types of grains: small and large 
Si-grains and small and large C-grains. The small Si-grain is simply referred to 
as small Si in this column. Each timescale is given in units of Gyr. The mark 
means that the listed physical process is not applicable. 
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Table 3 


The model 

PARAMETERS FOR 

RADIATION 

-DRIVEN 

DUST WIND 

IN EACH MODEL, 

Model name Cw,s,Si ^ 

C^w,l,Si 

Cw,s,C 

Cw.l.c 

comments 

SI 

0.01 

0.002 

0.05 

0.01 

fiducial SMC 

S2 

0.0 

0.0 

0.0 

0.0 

no dust wind 

S3 

0.01 

0.002 

0.01 

0.01 


S4 

0.01 

0.002 

0.03 

0.01 


S5 

0.01 

0.002 

0.1 

0.01 


S6 

- 

- 

- 

- 

no SB 

S7 

0.01 

0.002 

0.05 

0.01 

no SBl 

S8 

0.01 

0.002 

0.05 

0.01 

no SB2 

S9 

0.01 

0.002 

0.05 

0.01 

no SB3 

LI 

0.01 

0.002 

0.03 

0.01 

fiducial LMC 

L2 

0.0 

0.0 

0.0 

0.0 

no dust wind 

Ml 

- 

- 

- 

- 

fiducial MW 


^The first capital alphabets “S”, “L”, “M” indicate SMC, LMC, and 
MW models, respectively. The mark means no dust wind owing to 
no starburst events whereas the zero value (0.0) means no dust wind in 
spite of starburst. 

'’The dust removal coefficient for small and large Si-grains and small 
and large C-grains is represented by Cw.s.Si, C'w.i.Si) C’w.s.c, and Cw.i.c, 
respectively. 
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Fig. 1. — Time evolution of the total stellar mass (Mg), gas mass fraction (/g), and SFR for the MW (green solid), the LMC (red dotted), 
and the SMC (blue dashed). The SFR at each time step is the “normalized” one so that we can compare the time evolution of SFRs between 
the three galaxies more clearly. The normalized SFR is given in units of O.lMunit/^unit just for convenience, where Munit and tunit are mass 
and time units for a model, respectively, and they are initial gas mass (Mg) and 10® yr, respectively. 
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Fig. 2.— Extinction curves of the modeled SMC (black solid) at the selected 6 time steps, T = 8.4 Gyr (before the first starburst, “pre- 
SBl”), T = 9.8 Gyr (after the first starburst “post-SBl”), T = 11.5 Gyr (“post-SB2”), T = 12.5 Gyr ( “post-SB2”), T = 12.9 Gyr (during 
the third SB “SB3”), and T = 13.0 Gyr (“post-SB3”). For comparison, the observed extinction curves of the MW (green short-dashed), the 
LMC (red long-dashed), and the SMC (blue dot-dashed) are shown in each frame. 
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Fig. 3.— Time evolution of the total grain masses (normalized to the total mass of large Si grains at each time step for convenience) 
for small Si-grains green solid), large C-grains (Md(C,L); red dotted), and small C-grains (Md(C,S); blue dashed) for the SMC 

model SI with dust wind (upper) and S2 without dust wind (lower). Three thick yellow lines indicate the epochs of starbursts (SBl, SB2, 
and SB3). 
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Fig. 4.— Extinction curves for the eight SMC models. The left panel shows S2 with starburst (SB) yet no wind (green solid), S3 with the 
dust removal efficiency for small C-grains (Cw s c) being 0.01 (red dotted), S4 with Cw s C = 0.03 (blue short-dashed), S5 with Cw s C = 0-1 
(magenta long-dashed). The right panel shows S6 with no SB thus no wind (green solid), S7 with no SBl (yet with SB2 and SB3; red dotted), 
S8 with no SB2 (blue short-dashed), S9 with no SB3 (magenta long-dashed). The thick black solid line indicates the observed extinction 
curve of the SMC. 
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Fig. 5.— Same as Figure 2 but for the LMC model LI. The first starburst epoch (SBl) is different between this model LI and SI. 
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Fig. 6. 


Same as Figure 3 but for the LMC models with (LI) and without (L2) dust wind. 
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Fig. 7.— Extinction curves for the LMC models with dust wind (red solid) and without (blue dotted) at T = 13 Gyr (i.e, the present 
LMC). For comparison, the observed extinction curve is shown. 
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Fig. 8. — Same as Figure 2 but for the MW model Ml (without starburst and thus no dust wind). 
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Fig. 9. — Time evolution of dust grain masses (normalized by the mass of large Si-grain, Md(Si, L)) for small Si-grain (top), large C-grain 
(middle), and small C-grain (bottom), for the MW Ml (green solid), the LMC LI (red dotted), and the SMC Si models (blue dashed). 
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Fig. 10.— Dependence of the radiation pressure coefficient (Qpr) for silicate (solid) and graphite (dotted) with dust radii (a) of 0.01 
(top) and 0.1 /im (middle) on A and the spectral energy distribution (SED) of a galaxy with [Z/H]=—0.7 and age of 63 Myr (bottom). This 
SED with the adopted metallicity and age can mimic the SMC with a recent enhanced star formation. The frequency-averaged Qpr (QJr) 
silicate and graphite are indicated in the upper part of the top and middle panels. Clearly, Qpr is systematically larger in graphite tnan in 
silicate, which implies that graphite can be more strongly influenced by radiation-driven dust wind. 
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Fig. 11.— Time evolution of AMd^w(i)/^d(i) foi" Hi® SMC models with Cw =0.01 (green solid), Cw =0.03 (red dotted), and Cw =0.05 
(blue dashed) in the three starburst epochs, SBl (top), SB2 (middle), and SB3 (bottom). The same value for Cw (dust removal coefficient) 
is adopted for small and large silicate and carbonaceous grains (i.e., Cw,s,Si = C'w,i,Si = C'w,s,C = C'w,l,C = C'w) in these SMC models so that 
we can clearly demonstrate the physical meaning of Cw- The total mass of dust ejected from the SMC due to radiation-driven stellar wind 
(AMu w(^)) is estimated at each time step and AM^ w normalized to the total dust mass (Md(t)) at the time step is shown in this figure. 
Given the time step width of 10® yr in the present study, AM^ -^{t)/M^{t) = 0.01 means that the SMC can lose 1% of its dust through dust 
wind in 10® yr. 
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Fig. 12.— Same as Figure 2 but for the SMC model in which the strength of the first starburst (SBl) at T = 8.5 Gyr is twice as strong 
as SBl in the fiducial model SI. In this model, other model parameters are exactly the same as those of SI. 
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Fig. 13. — The time evolution of metallicities {Z normalized by the solar metallicity Zq) for the standard MW (Ml; green solid), LMC 
(LI; red dotted), and SMC (Si; blue dotted) models. 
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APPENDIX 

THE PHYSICAL MEANING OF Cw 

A key parameter for the time evolution of extinction curves of galaxies is Cw (dust removal coefficient) in the present 
study. Although an adopted value of Cw in a model can indicate whether dust can be removed more efficiently in the 
model than in other models with different Cw, the value of Cw (e.g., 0.01) itself does not represent the net rate of dust 
removal (e.g., O.OIMq yr“^) at each time step in the model. It is therefore important for the present study to relate Cw 
to the net rate of dust removal. Figure 11 shows the dust removal rate /Md{t)) at each time step in the three 

SMC models with different Cw. It is clear that (i) the dust removal rate is systematically higher for the model with larger 
Cw and (ii) it is higher in later starburst (i.e., higher in SB3 than in SBl). 

The dust removal rate of 0.01 at a time step during a starburst in a model means that 1% of dust at the time step can 
be removed from the galaxy within the time step width of the model (10® yr). This furthermore means that during the 
starburst with the duration of 10® yr, the galaxy can lose ^63% (=[1 — 0.99^°®] x 100%) of dust that it had before the 
starburst. The average AMd,w(t)/Md{t) is 0.0043, 0.013, and 0.022 for the three models with Cw = 0.01, 0.03, and 0.05, 
respectively. Therefore, AMd,w(t)/TId(t) (per Myr) can correspond roughly to 0.4Cw in the present study. It should be 
noted here that even for a constant Cw, AMd,w{t)/Md{t) of a galaxy can significantly change with time during a starburst 
owing to the rapid evolution of the total stellar luminosity of the galaxy during the starburst (See Figure II). The 
normalized dust removal rate (AMd,w(t)/fbfd(t)) might be useful for other theoretical models (e.g., one-zone models and 
numerical simulations with dust-related physical processes) that discuss the time evolution of dust in galaxies, particularly 
when they include radiation-driven dust wind. 

DEPENDENCE ON THE STRENGTH OF STARBURST IN THE SMC 

Although we consider that we have chosen a reasonable and realistic combination of Csbi, Csb 2 , and Cshs for the SMC, 
it could be possible that the present results could depend on these three parameters. A particularly interesting question 
is whether the SMC can have an extinction curve without the 2175 A bump yet with the steep rise in the FUV regime 
in the early evolution phase if there is a stronger starburst. We have accordingly investigated the SMC models in which 
only Csbi is different from the fiducial SMC model SI so that we can address this issue. Figure 12 describes the time 
evolution of the extinction curve in the model with Csbi = 0.1 (instead of 0.05 in the fiducial SMC model), which means 
that the first starburst is as strong as the third one in the model. 

Clearly, the extinction curves in the early phase (T < 9.8 Gyr) of this model are not so similar to the observed one of 
the present SMC, and this is true for the model even with Csbi = 0.2. These results accordingly suggest that an extinction 
curve similar to that observed in the present SMC can not be achieved in the present models at early evolutionary phases 
by changing the strength of SBl alone. These also imply that the SMC might have established its characteristic extinction 
curve quite recently. The final extinction curve in this model looks similar to the observed one for the SMC, but the 
degree of such similarity is higher in SI than in this model, which means that the present fiducial model SI chooses a 
more realistic parameter value than this model with Csbi = Csbs = 0.1 does. It is interesting that this model shows a 
more conspicuous 2175 A bump at T = 11.5 Gyr and 12.5 Gyr. It is confirmed that the observed SMG extinction curve 
can not be reproduced well in the SB2 phase of the SMG models with higher Csb 2 (=0.1). 

METALLICITY EVOLUTION 

The time evolution of metallicities (Z) of the MCs and the Galaxy has been already investigated in detail by Tsujimoto 
& Bekki (2009), Tsujimoto et al. (2010), and BT12. However, the models adopted in these studies are different from 
the present one in the sense that they did not include dust-related physical processes (e.g., dust growth). Therefore, it is 
instructive for the present study to present the results of metallicity evolution in the new models. Figure 13 shows that 
the global trends of the metallicity evolution in the present three galaxy models (SI, LI, and Ml) are very similar to those 
derived in our previous studies. However, the final metallicities at T = 13 Gyr are slightly smaller than the observed 
values of the galaxies: Z = 0.0046, 0.0058, and 0.013 for SI, LI, and Ml models, respectively. These Z slightly smaller 
than the observed values are due largely to the more efficient loss of dust and metals through SN and dust winds in the 
present new models. The models with less efficient SN and dust winds can show higher Z that are closer to the observed 
values. 



